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1. Introduction 



To accurately characterize imperfections of pre- 
cision coaxial air lines, skin effect and surface 
roughness need to be considered. Skin effect is now 
well documented [1] and conductor surface finish 
has been studied in detail by Rice [2] and Ament [3] 
through the use of Fourier series methods. While 
Karbowiak [4] points out that Fourier analysis re- 
veals useful knowledge of the spectral components 
which principally affect scattering parameters, it is 
also appropriate to examine local pointwise influ- 
ence along the axial (z) coordinate. In this connec- 
tion, Hill [5] developed perturbation expressions 
for the scattering parameters for a lossless circular 
air line. When the conductor surface exhibits trans- 
verse angular variation, Rouneliotes, Houssain and 
Fikioris report the effects of ellipticity and eccen- 
tricity on cutoff wave numbers [6]. 

The purpose of this paper is to develop numeri- 
cally accurate pointwise coaxial air -line scattering 
parameters that account for skin effect loss and 
conductor surface variations in the transverse an- 



gular and axial directions. Following Schelkunoff 
[7], Reiter [8], Solymar [9] and Gallawa [10], gener- 
alized telegraphist equations for the principal mode 
are derived in section 2 for a circular air line. 
Transformation to forward and backward wave 
differential equations enables general solutions for 
the scattering parameters in section 3. To allow for 
conductor surface measurements along the z-axis, 
cubic spline polynomials provide a starting point 
for establishing pointwise recursion formulas of 
forward and backward waves in section 4. In sec- 
tion 5, the Bergman's kernel technique is used to 
establish a conformal mapping for transforming 
noncircular conductors into equivalent circular 
conductors in correspondence to the principal 
mode. Computational results illustrating | S n \ ver- 
sus air-line length are given in section 6. An error 
analysis of the computational algorithms for the ac- 
curacy resolution of the measurement system is de- 
veloped in section 7. 
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2. Generalized Telegraphist Equations for 
the Principal Mode 

Consider the coaxial air line in figure 1. With 
inner radius a(z) and outer radius b(z) the field 
components of primary interest are radial electric 
field (E r ), angular magnetic field (Hg), and axial 
electric field (i? z ). We assume the fields E r , H g , and 
E z are composed of TEM and TM modes, and cou- 
pling of the modes is caused by skin effect with 
variations of the conductor surfaces. From Ap- 
pendix A boundary conditions for TM modes pos- 
sess the form 

E, = -Z,K{fo(zj}H 9 at r=b(z), (1.0) 

E, = +Z,K{4> a (?j}H, at/-=a(z), (1.1) 

where for instance, 

* {<Wz)> -l+tan [*,(*)]• ( } 

Appropriate Maxwell equations for determining 
transverse fields E r and He in the air dielectric re- 
gion of the air line are [1 1] 



dE r . TT dE z 

— = - Joi{ , Ha +— 



dHg . _ 

— = -jo>eE r . 



(1.3) 



(1.4) 



The parameters &), fi, and e are defined as radian 
frequency, permeability and permittivity, respec- 
tively, In addition the fields are assumed to vary 
with time according to the complex exponential 
function e*". 

To find the generalized telegraphist equations it 
is convenient to assume the fields possess orthogo- 
nal expansions in r and 6. In view of TEM and TM 
modes together with impedance boundary condi- 



-^H 




Figure 1. Coaxial air line. 



tions a set of orthogonal basis functions needs to be 
constructed from the Gram-Schmidt process. As- 
suming E r , H e , and E z possess continuous first and 
second derivatives implies their expansions are ab- 
solutely and uniformly convergent [12], In Ap- 
pendix B these properties are used to rearrange the 
expansions into the form 



E r (r,e*)= X Hi? CO e<£ (/-,e,z) 

(n,/>)=(0,0) 
(n,/>)=(0,0) 

+/S(z)A$(/-,0,z) 



(1.5) 



(1.6) 



where the superscripts (1) and (2) represent even 
and odd modes, respectively. We have 



.(i) 



(r,6,z) 



le™ {r,e,z)\ 



1 



cos kO 



~ N-^z/'^^X sin kO 



(1.7) 



In addition, 

a,A& (r,0,z)=a z Xa, < (r,0f) ; / = 1,2 (1.8) 
and N„ p (z) denotes the norm off mp (rz), that is, 



N np (z)= f Z {f mp {r,zYf mp (rf)y n &r, (1.9) 



where " stands for the complex conjugate. In par- 
ticular for the TEM mode 



e&(/v0 = 



1 1 



N 00 (z) = 



JVoo(z) r ' 

2. in^) 1 
a(z)J 



(1.10) 



Higher order modes are usual linear combina- 
tions of the first derivative Bessel functions /„' and 
Y' 

Following Reiter [8] by taking the inner product 
of eq (1.3) with the basis function e^ yields 

j ^ Si&dS = -jail J" zfle • {*&% XaJdS 
J^e^dS (1.11) 



*) 



+ 



&?) 



where S(z) denotes the cross sectional air dielec- 
tric region between the conductors. The left side 
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calls for differentiation of a variable surface inte- 
gral and the second member of the right side inte- 
grates by parts. 1 Hence, eq (1.11) evolves into the 
form 

SU) *) 

= § L tm{6 b (z)}E r e%ds 

- § L tm{<t> a (z)}E r e%ds -jafi f H e h e l 



dS 



S(z> 



+ § {E z [b (z),6j] -E 2 [a (z),6,z]}h%ds 



-J*2?,£s<£dS. 

S(z) 



(1.12) 



To express E z in terms of the constructed basis 
function (f„ p ), let the component of E z correspond- 
ing to k =0 (TM Qp modes) be expressed as 

£*=*£ b Qp (z) f /„» da. (1.13) 

p=\ J r 

From Maxwell's equations the (Op) mode relation 
between E. and E, is 






(1.14) 



where y /> and £ c , 0; , are the propagation constant and 
cutoff frequency numbers, respectively. Using eqs 
(1.5), (1.13), and (1.14) yields 



b Qp {z)=^N- p \z)V§{z). 

/Op 



(1.15) 



Now substituting basis function definitions eq 
(1.7) and calling for the principal mode yields 

dPoo " [ m f 0p jr,z) K „ 3 jUr^) 

= -jap mz)-2T Z s Kfaw fift^f biz). 
2^00 N 0p iz) 



1 Since S(z) is differentiable and the field components are con- 
tinuous, interchange of differentiation and integration is justi- 
fied. 



2,/W iV ,(z) Na&ifryotN^y 

fo P iu)du— jU.r,z)rdr 

a(z) J r °' 

-LO^rtc^rA <VM fooU> jz)^] h M y /op [b jz),z] ^, 
+2 W tan[<M*)] ^ (z) 6(z)X ^ (z) F^) 

-2. tan[cf, a (z)]^gfl « «t ^^gf 1 ^(z). 



^Voo(z) " v %4< N 0p iz) 



(1.16) 



To derive the companion generalized telegraphist 
equation from eq (1.4), the procedure is almost 
identical. Taking the inner product of eq (1.4) with 
h$ q from eq (1.8) yields 

j ^ h e % dS = -joe j a e E r • a 2 Xa,e ( ^ dS. 

5(2) 5(2) (1.17) 

For the left side 

i) ' =0 to 

^ (a, XM&l dS - § L tan {<j> A (z)} if„ /#£ dS 

+ § L tan {<f> fl (z)} He h 6 % dS. (1.18) 

Since integration by parts obtains the relation 

J e%± ~e% dS = - J e<^ e<J> dS, (1.19) 

5(2) 5(2) 

eq (1.17) takes on the form (setting k=q=0) 



d/po 
dz 



= -joe Vooiz)-2-7T J I Qp iz) f e<£^ e$ dS 



+ b'(z) f H e Affio dS -a'fc) #, #«, V&> dS. 

(1.20) 

Substituting eqs (1.6) and (1.10) into eq (1.20) gives 



/>=0 



W#£fJ£6?, 



J iVoo(z) 3z iV 0p (z) 



dr 



fl (2) 



yoi€Fffl(z) 









(2) ' 

(1.21) 
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Examining eqs (1.16) and (1.21) reveals that con- 
tinuous mode coupling occurs through the voltage 
and current transfer coefficients (left side), respec- 
tively, a phenomenon observed by Schelkunoff [7]. 
Skin effect coupling on the conductor surfaces was 
also reported by Schelkunoff and Gallawa [10]. 
When the air line is operated at frequencies appro- 
priate to the principal mode, all TM modes attenu- 
ate rapidly below their cutoff frequencies. 
Consequently, dominant coupling occurs between 
the forward and backward waves of the principal 
mode. 2 In this regard, eq (1.16) assumes the form 



d_Kx> 
dz 



r rbOD 

-2 J ■ 



1 



[ 



<&> 2ir In 






1 3 

1/2 - — 

r dz 



277-ln 



b ( z W /2 r dr V«fc)=—Jf>>p> Wz) 



a{z) 



z\k[^{z)] - 

1 b(z) In 



a(z) 



+*[*.(*)]- 



1 



a(z) In 



/ooCO 



+ 



b'{z) 



m. 

a(z) 
a\z) 



b(z)\n^\ 6(z)ln^. 
a(z) a(z) 



^Voo(z), (1.22) 



where the superscript, (1), has been dropped since 
only one mode is involved. Rearranging terms pro- 
duces the expression 



^ + {jap.+Z, k(z)}/oo=:Too W (z) Foo, 



where 



(1.23) 



<z)= mm + mg& i , ( , 24) 

b{z) a{z) j n £i£) 

a(z) 



m>(z)=\\ 



b'(z) fl'QQI 1 
b_ 
a(z) 



b(z) a(z)j ln 6(z)- 



(1.25) 



The equation for current proceeds similarly. Equa- 
tion (1.20) yields 



^ = -jaeV a >+T®Lz)I a > 



(1.26) 



2 Higher order mode influence on the TEM mode will be re- 
ported in a later issue of this journal. 



where the current transfer coefficient is defined as 



rffifcH 



b'(z) 



a'(z) 



Ku . . , b(z) . . , b(z) 
K b{z) In —■ a(z) In -fi- 1 
a(z) a(z) 



(1.27) 



3. Conversion of Generalized Telegraphist 
Equations to Forward and Backward 
Wave Equations 

Following Solymar [9] we define the amplitudes 
of the forward and backward waves A & and A ^ 
from the relations 

foo=«o \A oo+-^ oo J", 

■mx) = *o \A oo — A ooJ", 

where k =n/e, the wave impedance. Substituting 
eq (2.1) into eqs (1.4) and (1.20) produces the ex- 
pression 



I oo 



dz 



+ 



jfi+^^A&=T a ^)A& 



~^2£ A <» ~ 2"'ooCz)^( 



(2.2) 



dA 



oo 



dz 



jfi+ 



Z s k(z) 



2k 



Aoo=ToJz)A< 



Zk(z) + 1 + 

where /J=a)Vjxe, 

T<dz) =±{T®(z) + T®(z)}. 



(2.3) 



(2.4) 



In view of eqs (1.25) and (1.27) the last expression 
possesses the form 

T °°V-^b(z){b(z)-a(z)\- (2 - 5) 



a(z) 



For a lossless airline, voltage and current trans- 
fer coefficients assume the form, 



-TUz)=mz)=\— Jrc 
2. b{z) 



In 



b'(z) a'(z)} 



a{i) 



\b(z) a(z)\- 

(2.6) 
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Also, using Solymar's assumption that reflection of 
the principal mode does not affect forward propa- 
gation of the principal mode yields expressions 



&A\ 



dz 



and 



dAj 
dz 



+jfiA&=0, 



-jpAoo=- 



(2.7a) 



1 1 f fr'(z) _a'(z) 
2 ln W)\b(z) o(z). 
a(z) 



(2.7b) 



which agree with Hill's results [5]. 

Returning to eqs (2.2) and (2.3) and retaining 
Solymar's assumption above leaves the terms T m (z) 
A o>- Since coupling in this sense is meaningless, we 
drop the terms T w (z) A jjj, and obtain 3 



dAJo^l. nA Z s K(z) 
dz 

cL45b 
dz 



4< 



2k J 
Z s k(z) 



A&=0, 



2 Kg 



Am 



_\ Z s k(z 
1 2ko 



2) , T m (z) 



(2.8) 



(2.9) 



To incorporate appropriate boundary condi- 
tions, let the incident wave be A &(0)=Ao with per- 
fect termination at z=L, that is A^{L)=0. 

At this point the forward wave solution yields 



-U'W+ 



ZMQ. 



A&(z)= -Ae- u o w+ T^ wc> >0<z<L 



w , 



(2.10) 



and, at z=0, the reflected wave expression 

+ I^iy2ir wA adz (2 . n) 



show general forms which remain to be useful for 
using conductor radii measurements. 

From eqs. (2.6) and (2.7) the scattering parame- 
ters S n are S 2i are defined as follows: 



5ll =^and5 21 =^. 
Aq Aq 



(2.12) 



3 In the lossless air line 7oo(z)=0 in view of the eqs (2.5) and 
(2.6). 



4. Cubic Spline Fitting of Conductor 
Radius Measurements 

Underlying an accurate solution to A ob and A m 
are two critical items: (a) fitting conductor radii 
measurements with acceptable error bounds and 
(b) expansion of all known functions in a systematic 
manner to sufficient powers of z. 

To handle (a) consider cubic spline polynomials 
[13] for the inner (or outer) conductor measure- 
ments such that 



Q_i(z) = C ,k-i + ' • '-r-C^it-iZ 3 



(3.1) 



where & k _i(z) approximates a{z) or b{z), 
Zj t _i<z<z it such that k = \, • • •, N and Z N =L. It is 
desirable to transform the cubic spline eq (3.1) over 
the interval \z k _\,z k \ into the representation 4 



c*-i(0 = Qw-i+* • '+Q,it-i£ > 

in such a way that the condition, 

d 2 Q_.(D d 2 t 

d£ 2 -fei^-^*' 



(3.2a) 



(3.2b) 



holds at z=z k _i and z=Zjt where t,=z k —z k _ x . In 
addition we require 

Q_i(0) = & ,t-i> and 

Q;-i(zjt— Zjt_i)=& ,it-H r-Q.jt-iZjt (3.2c) 

such that C 0> jt_i represents the measurement of a(z) 
or b(z) at z=Zit_i. 

To implement (b), recall that E„ H e , and the sec- 
ond derivatives of E r and H e are assumed to be 
analytic functions in r, 0, and z. Hence, the expan- 

_1 can be 



sions 



of \b(z) ln^)- 1 and ja(z) ln^ 
I a(z)J I a(z)^ 



rearranged in powers of z. In Appendix C the fol- 
lowing expressions are derived over the interval 






(3.3) 



and 



4 A cursory inspection of calculated splines C* over the entire 
length of the line reveals that £bt is not necessarily equal to the 
radial measurement at z*_i. 
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2k ^ 



2 =2 8g,-.^, 



(3.4) 



where 0<£<z — z m _ l . 

To establish A m and A & at each point z n it is 
convenient to employ recursion relations. Inserting 
eq (3.3) into eq (2.6) and examining the interval 
Zi<,z<,z t yields 



^(0;z)=^3;z.)e-<C\V 1 ' , ' £ " dO 
where 
AMhzi)=A e-W l L'®>*«° 



(3.5) 



(3.6) 



is the forward wave emerging at z=z u 

For the interval z N _i<,zKz eq (3.5) immediately 
generalizes to the recursion relation 

AUO;z)=AUO?N-de- {f " N -\io S V»- ukd0 . (3.7) 

Proceeding to the backward wave A m by using 
eqs (3.3) and (3.4) in eq (2.7) for the interval 
z { ^z<,z 2 produces the relation 

A n(0?)=A mzd-A cWO^Ke-^'rV^ ""} 

■e- 2{/ oi * £Kd£} d7] (3.8) 

where the transformation t\=z—Z\ introduces the 
term e - ^ 2 '. Now eq (3.8) also generalizes to the 
recursion relation 

Aw(0^ N )=Am(0-^ N _i) 

-^^(O^Oe-^-'C^-'d 8 ^-^ 



3 m 1* +1 



(3.9) 



From Appendix C eq (3.9) assumes the solution 
A m(0;z n )=A m(P? n _0-A UO^-de-^-' 

•X vglv-i s k \-\yp+£ c® r _ l ]<? N -z N _X 

aTo L L ^o J J 

(3.10) 



where 



rtN 
Sk[a(zN— zn-i)]=\ ' 

•' Z N~\ 



e^z'dz. 



(3.11) 



5. Conductor Surface Variations in the 
Transverse Coordinates 

When the outer conductor is bored, circular 
cross sections are the exception rather than the 
rule. Most likely, an elliptical cross section evolves 
with some degree of rotation. Consequently, it is 
desirable to perform mechanical measurements of 
conductor radii in the transverse plane to charac- 
terize the deviation from circular cross sections. 
Since the principal mode is TEM in the transverse 
plane a direct conformal mapping of the measure- 
ment contour into an equivalent circular contour 
eliminates any difficulty of solving Laplace's equa- 
tion for an irregular boundary. If an equivalent cir- 
cular contour is found for each transverse 
measurement plane on the air line, a corresponding 
set of scattering parameters represents the original 
air line of measurement contours. 

The solution of Laplace's equation for a TEM 
mode with the inner conductor potential held to V 
and the outer conductor potential set at is 



, hi{r/b(z)} 

^- y °ln{a(z)/b(z)y 



(4.1) 



We initially state that Riemann's mapping theo- 
rem assures a mapping from the contour L to the 
unit circle and a particular expression for mapping 
evolves from the Bergman kernel expansion [14]. 
Thus, for a contour L centered at £ =8 e' 9 =0 the 
Bergman kernel is defined as 



5(o,o= Xi\(0) P V (D- 



(4.2) 



The Szego polynomials P V (Q are constructed to be 
orthogonal on the contour L and £ is a complex 
variable in the region bounded by L. Computing 
the inner (or outer) conductor radius requires the 
expression [15] 



c(z)= 



1 



Itt " 5(0,0) ! 



c(z)=a(z) or b(z), (4.3) 



where € defines the contour length of L at the 
point z. To see how eq (4.3) is constructed consider 
the differential line element on L, 

ds 2 =dx 2 +dy 2 , 

which in p, Q coordinates becomes 



ds 2 = 



p\0)+ 



-3£ 
30 



\d0 2 



(4.4) 
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Then the length of L is 

*-£*-£ Msir'"- (45) 

To find representations for the Szego polynomials, 
consider the following orthogonalization proce- 
dure. Let the matrix elements, 

hpq = ~k \ & 9ds= T f' ^ + ^> { P 2(0) 

+ [^]T e ' v " )9d *' (4 - 6) 

be defined for/>>0 and q>0. Note h pq =h qp . Then 
following Kantorivich and Krylov [15] compute 
the determinant 



D =\,D„ = 



"00 "io • • • "flO 
"o« "i« • • •"«« 
Hence, the Szego polynomial is defined as 



(4.7) 



p.®-- 



1 



[A-.A]' 



^00 h\0 ^n0 

Aoi. h u h„\ 

"0,«-l "1,»-1 • • •h'„,n-\ 

1 z • • • z" 



(4.8) 



A convenient property of the Szego coefficients 
for symmetric contours is found from eq (4.6). We 
have 



*""7".r {'-H-irwVff) 



+ [*],} e<— d». 

For off diagonal elements 

h M = 0;p+q odd, p=£q. 

Hence, any asymmetry in the contour L is expected 
to be noticeable through the off-diagonal elements 
h„. 

To find the equivalent circular conductor radius, 
integrals in eqs (4.5) and (4.6) need to be deter- 
mined from measurements of p and on the con- 
tour L . Let the following cubic spline be defined 



p,_,(0)=XPa-.e' 



(4.11a) 



such that Pk-i(0) approximates a(Q,z) or b(0^) 
over the interval 0*_i<0<0* for k = \,—, N, i.e., 

t7Ar = 277\ 

Following the same procedures as in the trans- 
formation from z to £ in eqs (3.2a) to (3.2c) enables 
the cubic spline, 



such that 



-j-J L P„(QP m (Odr > =S„ 



(4.9a) 



Carrying out the above procedures yields the con- 
ductor radii accurate to third order, 



f^ - ' 
C_ 27T 


i "10 #10 

AA 


(^io^i — ^n^2o)(h 10 h2i — huh 2 o) 1 


A AAA r 


where 


D\ = h n —ho\h w , 


and 


D 2 =(h n h 22 — h a h 2i )— hi (hoih 2 2— h 02 h 2i ) 


+h 20 ( 


^0lAl2 — Will)- 



(4.9b) 



(4.10) 



P;fc_l(<J>) = 2(P/,*-l)<J>', 

(=0 



(4.11b) 



to be constructed where &=0 — Q k _ x and po,*-i 
equals the measurement of a(0^) or b(0jz) at 
= 0t_i. 

Equations (4.5) and (4.11) yield an expression for 
length using the binomial expansion: 



<f 



=f"[p(0)+\ P -Xe) 



3£ 

dd 



d0. 



(4.12) 



Appendix D, taking into account the spline coeffi- 
cients, produces 

N-i re„+l-«„ f 3 7 1 

/ »=I 2 p ^ e *+5 2 p S- e * \ de > 

(1=0 J U=0 /t=0 > 

(4.13) 



«=o U=o a. -r i 

+ 2 2/^» *+i ) 



Expressions for the coefficients A w are developed 
also in Appendix D. 
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6. Computation Results 

The amplitude of <S U has been computed from 
eqs (2.12) and (3.9) for a 7 mm air line approxi- 
mately 15.6 cm in length using the frequencies 6, 
12, and 18 GHz. In addition the number of conduc- 
tor dimensional measurements in three sections of 
air line with variable spacing is shown in figure 2 
and the results are illustrated in figures 2-5. Mean 
and standard deviation values of the conductor di- 
mensional measurements are as follows: 



Inner radius Outer radius 



Mean (meters) 

Standard deviation (meters) 



0.1521 XlO" 2 
0.2481 XlO" 6 



0.3500X10" 2 
0.6836X10- 6 



Figure 2 illustrates that conductor radius measure- 
ments near either end are more volatile — particu- 
larly the outer conductor. Figures 3, 4, and 5 reveal 



that changes in conductor radii in the z-axis 
provide the dominant contribution to \S n \ while 
skin effect loss amplifies the in-phase and out-of- 
phase behavior of the lossless air line (as shown in 
fig. 2). In addition, skin effect loss affects the most 
significant digit of |5n| even for short lengths of 
line. On comparison of figures 3, 4, and 5 with 
Hill's results [5], the most noticeable feature is the 
overall difference in magnitudes of Su, which 
evolves from a uniform inner conductor model and 
lossless boundary conditions in Hill's work. Con- 
formal mapping effects from elliptical measure- 
ment contours do not affect 5 M and S 2 i unless the 
eccentricity is greater than 5 XlO -6 meters. How- 
ever, if the inner conductor has an eccentric posi- 
tion with respect to the outer conductor, 
conformal mapping by Bergman's kernel reveals 
scattering parameters Su and S 2 i are noticeably af- 
fected. 5 
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Figure 2. \Sn\ vs variable length of 7 mm lossless air line. 



s For additional computational results see Holt [17]. 
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Figure 3. \Su | vs variable length of 7 mm air line with skin effect loss at 6 GHz. 



7. Error Analysis 
7.1 Error Sources 

Error sources that contribute to scattering 
parameters evolve from (a) spline interpolation 
with respect to z, (b) spline interpolation with re- 
spect to 0, (c) conforfflal mapping using the 
Bergman kernel, and (d) expressions for backward 
and forward waves. 

To examine (a) and (b) consider the error bound 
from cubic spline interpolation theory [13], 

\f(x)-SA k (x) | <| [AI }i(f{x); | A, | )], (6.1) 

where SA k (x) defines a cubic spline, /i, signifies the 
modulus of continuity, 6 /" denotes the second 
derivative of the function /, and A* stands for the 
mesh size between arguments x k . An approxima- 



tion to the modulus of continuity is 



,*[/•"(*); | A* | ]=max 



/[jCfc-l,JCfc] 

h 



h 



2h 



(6.2) 



where f[ ] denotes the first order divided differ- 
ence of/. 

For a 7 mm air line with a mesh size (A*) equal 
to 1 mm, a value of p, from observations of conduc- 
tor radii measurements as functions of z indicate 
ii,=0.1 is reasonable, and in the angular direction 
/i,<0.01. 7 Therefore, the total error from spline in- 
terpolation is (considering the errors as additive) 

Error totaI = Error 2 + Error«<2.8 X 10" 7 m. (6.3) 



1 See Davis [16] for a suitable definition. 



7 Since mechanical measurements for roundness are not avail- 
able, an ellipse of 38.1 X 10 -6 m was selected for determining the 
modulus of continuity. 
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Figure 4. |5n | vs variable length of 7 mm air line with skin effect loss at 12 GHz. 



The error source (c) from conformal mapping an 
elliptical measurement contour of eccentricity 
equal to 38.1 X 10 -6 m (150 ju, in) is illustrated be- 
low. 

First term of outer conductor equivalent 
(mapped radius) 3.4989 mm. 

Second term 2.32X 10 -12 mm. 

Third term -3.35X10 -5 mm. 

Keeping in mind that ellipses are symmetrical 
with respect to the origin reveals that even terms 
are effectively zero (on the order of 10 -12 in view 
of machine precision). Since the convergence 
above is very strong, the fifth term is likely on the 
order of 10 -9 mm. The error; from source (d) de- 
pends on the number of expansion terms represent- 
ing the functions l/a(z), \/b(z), and ln{b(z)/a(z)}. 



At this point accuracy considerations of the for- 
ward and backward wave eqs (3.7) and (3.10) in 
correspondence to the measurement system are in 
order. For instance, to examine A^, consider the 
measurements of a(z) and b(z) at z = and z=z { . 
From eqs (3.7) and (3.10) we have the total differ- 
entials 



AA&(PZi)=A 









a b 



and 
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Figure 5. \Su | vs variable length of 7 mm air line with skin effect loss at 18 GHz. 



AAUO*0= -A [§ o UUaM^ [£j - *]) 



a 



^-2W<>a,bo.& 



+- 



In 



a, 4,11 e -2fft^,^)z 1 _ 1 
a 2 o 6 2 oJi" 2W(fl ,b^ s ) A/ "° 



+^ ^-{C/i,.(ao,6o) 



1 



In 



bo 



1-2]K**» 



Uljfioybo) 



1 



ao+VlJl^'W'o^) 
1 



W(.a*b Z,)' r W 2 (a ,!H>Z,y 



Wtao,Ao.Z>l 



A/b 
(6.5) 



where A/- is the measurement precision common to 
a and b . In the expressions above 



U^ ik {a ,b )——j+-rj 

"0 °o 



In* 



(6.6) 



and 



JJWo.Z,) =|^ C/ u ( ao ,6o) + {fj - ^} (6.7) 

where ai and bi are the coefficients from differenti- 
ating the cubic spline representations of a(z) and 
b(z), respectively. 

In computing AAqo and AA^ for a 7 mm 
air line let the measurement precision be 
Az=Ar =2.8xl0 -7 m for a frequency range of 
1-18 GHz. Using the measurements of a(z) and 
b(z) in figure 3, we select the maximum divided 
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difference magnitudes for a t and b\ to obtain the 
results 



I A^S,(0;z,)| =3.6X10" 5 at 1 GHz, 

= l.lXlO" 4 at 18 GHz, 
|A^5b(0^i)| =4X10" 8 at 1 GHz, 

= 1.5XlO" 7 at 18 GHz. 



to a conductor surface measurement resolution of 
2.8 XlO -7 m over frequencies appropriate to the 
principal mode for 7 mm air lines. 



(6.8) 9 # Acknowledgments 



(6.9) 



To examine the total uncertainty in A &> over N 
measurements of a(z k ) and b(z k ) for k = \,—,N let 
AA M be represented as 

±A &(0?n) = 2 f<- a (Zk)MZn)}±A U*n- I A), 

(6.10) 

where /„ is on the order of Aw(z n -\\z n ). Since the 
computation of LA^(z n _^ n ) proceeds in the same 
way as AA^(0]z{), in eq (6.5), the total uncertainty 
over all measurement positions is found when the 
individual uncertainties in z and r at each measure- 
ment position are known. 

7.2 Efficiency Improvements in Cubic Spline 
Approximation 

While fitting the surfaces of coaxial air line ge- 
ometries with products of cubic splines over the 
variables 6 and z successfully meets error bounds 
consistent with measurement precision, significant 
reductions in the number of measurements yields 
equivalent error bounds with Gordon's successive 
decomposition spline [18]. The number of measure- 
ments required for successive decomposition 
splines in comparison to usual spline products is 
generally less than 50 percent. 



8. Summary 

Generalized telegraphist equations for the coax- 
ial air line have been derived under two assump- 
tions: (a) skin effect losses are present, and (b) 
conductor surface variations occur in the axial and 
transverse coordinates. Product cubic spline ex- 
pressions to accurately fit conductor surface mea- 
surements were employed to arrive at pointwise 
scattering parameter expressions. Error bounds 
from eqs (6.8) and (6.9) reveal at least four signifi- 
cant figures can be obtained to characterize the 
scattering parameters S 2) and S u in correspondence 
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Appendix A. Inner and Outer Conductor 
Boundary Conditions 

In reference to figure 6 the boundary condition is 

E v =zJ v (A.l) 

where E v is the electric field component, J v is the 
current density component in the coordinate v, and 
Z, is the radial impedance at r=b(z). 
Also from figures 1 and 6 



a v = | a v | cos{<f>i(z)}a z +\a v \ sin{<f> 6 (z)}(-n), 
a v =cos{<f>i(z)}a z — sin{<f>4(z)}n, 



(A.2) 



where n is the unit outward normal vector. From 
the left side of eqs (A.l) and (A.2) 

E v a v = | E v | cos{<f> 6 (z)}a z —\E V \ sin{<f> 6 (z)}n, 

E v a v =E z a z +E r a r . 




Unit vector an is 
tangent to inner 
surface of outer 
conductor at 
(b (z), z) 



Detail of coordinate system 



Proceeding with the right side of eq (A. 1), yields 
J v a v =J z a z — J„n=a z • nXH— J„n. (A.4) 

To find an expression for J n note that 
J n =n- J„=n-n t xH, (A.5) 

J n =n • n^X {H T a T +H n n+H z a z }, 
J„=n • n ()J x{-fl r Tar+-fl r „{-sin{<f> i (z)}a v 

+cos{<f>4(z)}n^}. (A6) 

+# z {cos{<f> i (z)}a v +sin{<f> i (z)}n (t J, 
/„ = — sin{<f>i(z)}H T . 

Consolidating eqs (A.3), (A.4), and (A.6) yields 



E z = —Z s 



l+sin{<j>f,(z)> 



l+tan{<fo,(z)}. 



Ha 



(A.7) 



where the relation E z tan<f>(z)=i? r has been em- 
ployed. 

The inner conductor boundary condition rela- 
tion is similarly derived as follows: We have 
E v —Z s J v and the left hand side is 



E V =E Z +E r . 
The right side becomes 

av*/v = ^ir z !&/*'/■• 

Now proceeding as before yields 
J r =a r -J v =a r -n^xH, 

J r ^ 3 r • tiq B. vt 

J r =sin{^> a (z)}H g . 



(A.8) 



(A.9) 



(A. 10) 



(A.ll) 



Substituting i? r =tan{<J> a (z)} E r into eq (A.8) and 
using eq (A. 1 1) in (A.9) yields 



l+sin{4>„(z)} 
£ *- Z M+tan{«Mz)}^" 



(A. 12) 



Figure 6. Outer conductor coordinate systems for determining 
boundary conditions. 
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Appendix B. Transformation of Basis 
Functions for the Field Components 



Fg$Cz)=i>fc/.QCO^CO- 



/=0 



The transverse fields E r and H„ are represented Proceeding to the rath term yields the expression 



as 



E r (r,Q,z)=2 f <t> v (rjy{P$<?) cos(« 0) 

+ ft«(z)sin(rt0)} (B.l) 

H e (r,e,z)=2y r *„(r,z)GgOO cos(n9) 



HS(z)= 2 a Ay>ra (z) K<?) (B.5) 

for the coefficient of the basis function f km (r,z)/ 
N km (z). Now eq (B.4) from eq (B.5) becomes 



/*2tt 

E r (r,e,. 

J 



,z) cos&fldfl 



+?$(z)sin(«0)} 



(B.2) 



where <1>„ P denotes the radial component of the nor- 
malized potential function at a point, z. The symbol 
* over V„ p (or I„ p ) signifies an amplitude corre- 
sponding to the nonorthogonal function — <b np and 

the superscripts (1) [and (2)] correspond to even 
and odd modes. We assume that E r and H e have 
first and second derivatives in r and 6 which are 
continuous and bounded. Hence, the Fourier series 
eqs (B.l) and (B.2) are absolutely and uniformly 
convergent; that is the series can be rearranged. 

Suppose the kth even mode of eqs (B.l) and 
(B.2) is considered and let the Gram-Schmidt or- 
thogonalization be represented as 



p 
I 



- 4>„ p (rj)= £ a^-us (z)f k , (rj)/N k ,(z) (B.3) 



where f k e(rj)/N k f(z) is the k <fth orthogonal basis 
function. Now substituting eq (B.3), multiplying eq 
(B. 1) by cos(k0), and integrating over 9 yields 



f 

J o 



E r (r,02)cos(k8)&e 



= i P*o(z) a*,o,o(z)/Ao(r ) z)/WAo(z)} 
^^W^/^ + a k ^/^} + --- 



+ ^(z)K m .o(z)^g 



+ -+a k , m , m (.z)-^^}+..Mfc). (B.4) 

By induction note that the basis function J *o(r,z)/ 
Nkd(.z), occurs in every term above. Hence, we set 



=v(*) X VWUirtyN^z) (B.6) 

p=0 

for the kth even component of E r . The odd com- 
ponents of E r are obtained similarly. Hence, sum- 
ming all components of E r obtains the desired 
result 

E r (r,0*)= v(n) £^S W z > cos( - n ® 

n,p ly "P\ Z ) 

+ Oz)sin(«0)} (B.7) 

where v(«)=2tt if n =0 
=7rif n=£0 

The expansion for H e is obtained in the same way, 
i.e., 

H e (r,e,z)=v(n) £^S G#(z) cos(«0) 



+/£>(z)sin(«0)}. 



Appendix C 



(B.8) 



Beginning with the spline expressions eq (3.1), 
consider the reciprocal 



1 1 



1 



b(z) b . bi h, 

i+, z + »+, z 
b bo 



We define the coefficients 
Bf ) =b 1 /bo,-,Bf>=b 3 /b . 
Since the geometric series for small z, 



(CI) 



(C.2) 



0o b 



<a <<l 
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converges uniformly and absolutely, eq (C.l) can 
be rearranged in powers of z. We have 



1 1 



1 



b(z)-b \+B<?z + ~+Bfz>- 
Using the geometric expansion yields 

—=j-{\-BVz B^ + itBffz 2 



b(z) b 

+2B\°^z i +...}-[Bffz i —.}, 



(C.3a) 



and rearranging the above series into successive 
powers of z produces 



1 -=f +£>$?* +MV+MV- 



b(z) b 
where 

D% B\°\ D®= -Bf+[B\y 

D%= -Bf+2 BV B^-{B?n\ 



(C.3b) 



(C.4) 



To obtain an expression for ln{b(z)/a(z)} consider 
the expression 



In 



b(z) 



=]n^+\n{l+B\ 0) z+-+Bfz 3 } 



-ln{l+A\°>z+~+Afz 3 } 



(C.5) 



Taking the reciprocal of eq (C.5) and using the ge- 
ometric expansion yields 



In 



\bjz} 



[a(z)\ 



h SM i -{ h £H toI1+J * 



+..]+ln[l+A\Vz+..] 



+ 



ln^ 



2 {ln 2 (l +B?>z + ~) 



-21n[l+5\ 0) z + ..] ]n[l+A\ 0) z + ~] 
+ln 2 [l+^rz + •■]}. 



(C.6) 



From the absolutely and uniformly convergent ex- 
pansion 



. r. i O) 2 CO 3 

ln{l + 0)} = 0)-y+y + - 



for |a>|<fl <<l 

the expansion eq (C.6) is rearranged in powers of z 
to obtain 

w-i*" <c7) 

where the first four coefficients are defined as 

"oj 

! )" 2 {G|P-Gi!>} 



[ ln^ 



a 






r a i- 2 

a 



+ 



In 



b 



a \ 



H$= - 



(C.8) 



lnM" 2 {G|P-G<P} 
+ {ln |}" 3 .{Gg>-2Gff G<i>-2Gj£> G<i>+G<?} 

- {in ^}"*{Gg>-3Gg> G#+3Gfi> G£>+G#} 

and the G coefficients for the outer conductor cue 

G#={flr}*;* = 1,2,3. 

GJP=M°>-^|^, 



Gf?=2Gg> G?>=lBf 



M 0) - 



w 



(C.9) 



G^-B^-B?*^ 



[5H 3 



The G coefficients for the inner conductor are 
computed by replacing B k (0) with A k (0). To obtain 
the expansion for k in eq (1.15) an expression for 
K[<f>(z)] needs to be found. From eqs (3.1) and (1.1) 



r , ,, l + sin[4> t (z)] 



+tan[<f> A (z)] 



(CIO) 



l+6'(z) 
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Using the binomial and geometric expansions, in and ju.(«)=0 for n <3, 



the last expression, yields the following expression 
in terms of spline coefficients 

K[^ b {z)\ = \-\b\-6b,b 2 z. 



We now have from eqs (C.3), (C.7), and (CIO) the 
expression 



K{4»(?)} 



£>(z)ln 



b(z) 
a(z) 



-H 



b\ 



—6bib 2 z 



£JD^"2#<!?z™=2Cfi>z* (C.ll) 
oo o 

where C ( ^ is the convolution of coefficients, i.e., 

cjp=(i -\ bijiDm+Hmy 

-6b,b 2 D^H^ (C.12) 

CS>=fl -| b{) ZD±°l k M" -6 b,b 2 ^ D\* k m> 

V L /*=0 *=0 

Since the inner conductor expansion proceeds in 
the same way, the total expression is 



K (z) = ^Cfz k (C.13) 

k=0 

where 

cV=c$+c®. (c.14) 

The transfer coefficient eq (2.4) has the expansion 



T M (z) = [b, + 2b 2 z + 3b 3 z i i [£ £>#z*] [£ ^M 
-[a l +2a 2 z+3a 3 z 2 ]jl D ® z " X ^ z *' 



*=0 



*=0 



r M (z) = £ C?>z 



?),* 



(CI 6) 



where 



n m 



_(„ + 1 _ m ) a „ +1 _ m jDm_ tJ ffj t 1) } 



= 1 for «=3. 
From eqs (2.3), (C.14), and (C.17) set 

8(f = M C?+ 2 ;k> °- 



(C.17) 

(C.18) 
(C19) 
(C.20) 



The coefficients vk%_i in eq (3.14) are defined by 
selecting the coefficients of rf and convoluting the 
series product. Let 

vSK.-i=i 

v&-i=0 



va-^-^Sg^.+^Si".,} 2 



vgL-^fSft,-,^., 



vJ»,-.=2{ i % ji +|8i'jL_i8g),-,' 



J,m-1 



vgJ-^IJs^-i} 2 -^^,' 



Proceeding to vjt^-i 

where 

/Xi(fc)=&if 0<&<3 
= 3 if 3<A:<8 

and 



(C.21) 



(C.22) 



(C.23) 



v@L-i= 2 ^ 2 > Am _,v<J> m _ 1 if8<fc<ll. (C.24) 

€=k-Z 
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Appendix D 

Expression of the contour length <? in terms of 
the spline coefficients is initiated from the defini- 
tion 



/ =J 2,r {s(0)+r 2 (0) 



dd. 



(D.l) 



Using the spline coefficients for the interval [0„, 
Q„+i] yields the representation 



Pn(v)=^PilnV k ,0<V<Gn +i 



(D.2) 



k=0 



at n =0, P n (0)=P\% n =the measurement of £ at 
0=0„. Then eq (D.l), in view of eq (D.2), yields 

n-i /-e„+i-e n r 3 



„=0 J lfc=0 



+^(Hl)4i/ 



fc=0 



]d0.(D.3) 



Let P&' " be the convolution of the product 
*& + »«- l5 and {(* + 1) P$ +hn } 2 above. Then eq (D.5) 
yields 



1 



3(p+q+l) 



Ko=-Tj J ^- qy> " £ ^& ?+1) {s,[(0 B+1 -0„) 

■* fc=0 



■■Cs*[(^ + i-e»Kp-?V]-5*(0)} 

where 5 k is defined in eq (3.11). 



(D.7) 



The denominator can be expanded into a geometric 
series such that the second term above yields 



i ?st,0*=f i (*+i) m + u„e k \ 2 \4r 

k=0 U=0 J l-MA* 



(D.4) 



and eq (D.3) on substituting eq (D.4) yields the 
form eq (4.13). For other values ofp and q we have 

^ o «=o Jo IS 

k=0 i U=0 



+ 1/2 



r3(p+ ? +l) 

X 

L fc=0 

.e*'-« )e d© 
where 



(P+1) f 3 l/>+? 

X p&e k =\2Pi!lo k \ ■ 

k=Q U=0 J 






(D.5) 



(D.6) 
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